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ABSTRACT 

We consider solitons in a nonlinear optical fiber with a single polarization in a region of 
parameters where it carries exactly two distinct modes, viz., the fundamental one and the 
first-order helical mode. From the viewpoint of applications to dense- WDM communica- 
tion systems, this opens way to double the number of channels carried by the fiber. Aside 
from that, experimental observation of helical (spinning) solitons (that can be launched 
and detected, using helicity-generating phase masks) and collisions between them and 
with fundamental solitons in (ordinary or hollow) optical fibers is an issue of fundamental 
interest, especially because it has been very recently found that spatiotemporal spinning 
solitons in bulk optical media with various nonlinearities are unstable. We introduce a 
system of coupled nonlinear Schrodinger equations for fundamental and helical modes, 
computing nonstandard values of the cross-phase- modulation coupling constants in it, 
and investigate, analytically and numerically, results of "complete" and "incomplete" 
collisions between solitons carried by the two modes. We conclude that the collision- 
induced crosstalk is partly attenuated in comparison with the usual WDM system, which 
sometimesamay be crucially important, preventing merger of the colliding solitons into a 
breather. The interaction between the two modes is found to be additionally strongly sup- 
pressed in comparison with that in the WDM system in the case when a dispersion-shifted 
or dispersion-compensated fiber is used. 
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1 Introduction 

A commonly adopted approach to the description of nonlinear light propagation in optical 
fibers is based on separation of the transverse modal structure, that may be described 
in the linear approximation, and slow longitudinal and temporal evolution of the signal's 
envelope, which is essentially affected by the temporal dispersion and Kerr nonlinearity. 
The analysis following this way ends up with the derivation of the nonlinear Schrodinger 
(NLS) equation for the envelope |], ||. 

Usually, both experimental and theoretical studies of the soliton propagation are con- 
fined to the case when parameters of the fiber admit a single (fundamental) transverse 
mode, because in a multimode fiber an initial pulse excites different modes in an uncon- 
trollable fashion. However, using well-known data for the fibers of the simplest step-index 
type [Q, it is easy to find that a situation with exactly two modes takes place when the 
standard waveguide parameter 

V = kpy/nl - n% (1) 

where k, p and n CO)C i are, respectively, the carrier wave's propagation constant, core radius, 
and the refractive index in the core and cladding, takes values 

2.405 < V < 3.832 . (2) 

For instance, in the case of the standard carrier wavenumber A = 1.54 //m admitting the 
soliton propagation in optical fibers, and taking the usual value n co — n c \ = 0.01, the 
interval (0) implies 3 //m < p < 4.75 fim, i.e., quite realistic values of the core's radius. 

Inside the interval (|2|), the fiber carries a fundamental mode (FM) and the first helical 
mode (HM). The transverse structure of the latter one is described by the expressions 
J\(Ur) exp(±z0) in the core and K\(Wr) exp(±z#) in the cladding, where J\ and Ki 
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are the standard cylindrical functions, U and W are the associate waveguide parameters 
defined in the usual way [[J, and (r,9) are the polar coordinates in the fiber's cross section. 
In this case, one is actually dealing with a set of three coexisting modes, as there are two 
degenerate HMs with the helicities S = ±1. Note also that, because the physical fields 
are proportional to the real part of the complex expressions, the presence of the multiplier 
exp(±z6 l ) means that HM solitons are spinning in the course of the propagation along the 
fiber. 

As neither higher-order radial (nonhelical) mode, nor any HM with a helicity > 1 
exists in the interval the two-mode situation is controllable: a light pulse with zero 
helicity can excite solely FM, while excitation of HM is possible by a pulse that carries 
the necessary helicity. A light beam can be lent helicity, passing it through a specially 
designed phase mask, which is quite feasible in a real experiment, see, e.g., Ref. ||. Due 
to their distinct topological nature, FM and HM do not linearly mix, provided that the 
fiber remains circular. As is well known, it is easy to fabricate a long silica fiber whose 
deviation from circularity is negligible. Fiber's bending will induce no linear mixing 
either, providing that the bending radius is much larger than the wavelength. In this 
work, however, we do not discuss exact limitations on the deviation of the fiber from the 
circularity and similar details. Instead, we focus on principal issues, such as collisions 
between solitons carried by FM and HM. 

It is necessary to distinguish between the fundamental and helical solitons at the 
receiver end of the fiber. In an experiment with a single or a few copropagating carrier 
frequencies, this is quite simple, as the fundamental and helical modes have an appreciable 
difference in their propagation constants (see below), thus the two types of the solitons 
can be discriminated between by means of a simple wavelength filter. Besides that, there 
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is a possibility to create a "helicity filter" , which would also work in the case of a multi- 
channel WDM system. Indeed, if it is known that two species of the solitons in the fiber 
have the helicities S = and S = +1, at the receiver end the incoming signal can be 
passed through a phase mask that adds extra helicity AS = —2. Then, the arriving FM 
and HM solitons will change their helicities to —2 and —1, respectively, so that only the 
latter one will survive, as in the selected parametric region the modes with S = ±2 do 
not propagate in the fiber. On the other hand, passing the incoming signal through the 
phase mask adding AS = +1 will transform the former S = soliton into a propagating 
S = +1 pulse, while the former S — +1 soliton will have 5* = +2, hence it will not be 
able to propagate. 

Thus, launching solitons independently in each of the two modes, one can implement a 
two-channel system inside the core. Note that standard elements of the fiber communica- 
tion systems, such as amplifiers and guiding filters 0, will act in essentially the same way 
on the solitons in both modes (although the gain coefficients of an Er-doped or Raman 
fiber amplifier may differ for the two modes, depending, e.g., on the density distribution 
of the doping atoms in the fiber's transverse plane). Moreover, if one starts with a WDM 
(wavelength-division-multiplexed) multichannel system already implemented in the fiber, 
one can double the number of the channels by means of this two-mode scheme (we demon- 
strate below that it is not really possible to triple the number of the channels, using two 
HMs with the opposite helicities). The feasibility of such a "mode-division" channel dou- 
bling may be quite important, as it has been demonstrated that doubling by means of 
the polarization-division multiplexing is incompatible with WDM |J. Indeed, while the 
polarization of the soliton can be easily changed by various imperfections of the system, 
the mode's helicity is expected to be robust, as it is a topological invariant. Note that we 
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do not consider the polarization structure of the modes, assuming that either they belong 
to one polarization, or (more realistically for the applications) the polarization can be 
effectively averaged out. Thus, our helical mode has nothing in common with the circular 
polarization. 

Due to the Kerr nonlinearity, the linearly orthogonal solitons borne by the two modes 
interact via the cross-phase modulation (XPM). The main technical objective of this work 
is to study XPM-induced effects of collisions between the solitons. As is well known, the 
collisional crosstalk is the most fundamental problem in the soliton-based multichannel 
communication systems, see. e.g., Refs. f| - §§. 

It should be stressed that, while the application of the proposed mode-division dou- 
bling to WDM soliton communication systems is not straightforward, as some technical 
problems remain to be resolved, experimental observation of narrow subpico second he- 
lical (spinning) solitons and their collisions between themselves and with fundamental 
solitons in relatively short optical fibers is a problem of fundamental physical interest by 
itself. A combination of the above-mentioned helicity-generating phase masks with the 
well-developed experimental technique admitting, e.g., direct observation of the polariza- 
tion structure of subpicosecond solitons in short fibers || should make the observation 
of the helical solitons and their collisions quite feasible. An additional interest to the 
latter problem is lent by the fact that a similar object in a bulk optical medium, viz., a 
spinning light bullet (spatiotemporal soliton), has been very recently found to be unstable 
in models with various nonlinearities [IIJ. Therefore, the optical fiber may be a unique 



medium in which the existence of a stable spinning temporal soliton may be possible. 

It is also noteworthy that the helical soliton, whose local intensity vanishes at the 
central point of the fiber's cross section, may be a natural object to exist in hollow 
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nonlinear optical fibers, which have recently attracted a lot of attention (and where, 



incidentally, very narrow solitary pulses are quite possible), see, e.g., the works [11] and 
references therein. An interesting issue is a possibility to select parameters of the hollow 
fiber so that it would support solely a helical mode, which would then play the role of the 
fundamental one. 

The paper is organized as follows. In section 2, a system of coupled NLS equations 
for the three modes, FM with S = and two HMs with S = ±1, is obtained. It has 
nonstandard values of the XPM coupling constants. In section 3, collisions between 
solitons carried by the fundamental and helical modes is studied analytically, by means 
of the perturbation theory. The perturbative treatment applies to the case when the 
colliding solitons pass through each other quickly enough. In section 4, an example of 
direct numerical simulations of the collision between the FM and HM solitons, illustrating 
a difference from the collision between FM solitons in the usual WDM system, is displayed. 
The difference may be crucial: the colliding solitons merge into a breather in the usual 
system, but survive the collision when they belong to the different modes. The paper is 
concluded by section 5. 

2 The Model 

A normalized system of coupled nonlinear Schrodinger (NLS) equations for the interacting 
modes can be derived by means of a standard asymptotic procedure |J, 

i(u ) z + ic (uq) t + k u - ]-(3 Q (u Q ) TT + (|u | 2 + 27 |u+| 2 + 27 |m_| 2 ) u = 0, (3) 



i(u + ) z + id (u+) T + k x u + - -/3i(tt + ) rT + (|w+| 2 + 2|w_| 2 + 27i|w | 2 ) u + = 0, (4) 
i(u-) z + id (u-) T + hu- - ^-(3i(uJ) TT + (Jm_| 2 + 2|m + | 2 + 27i|m | 2 ) u~ = (5) 
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(in the case of extremely narrow solitons, well-known higher-order terms [p], @] should 
be added to the system). We here consider the most general case, when two HMs with 
the helicities ±1, represented by the envelopes u±, interact with the zero-helicity FM ii ; 
/?o and Pi are the corresponding mode-dependent dispersion coefficients (see below), co,i 
and fco,i are the group-velocity and propagation-constant shifts of the two modes (these 
characteristics are also mode- dependent ||), and the effective XPM coefficient 7 and 71 
are given by the properly normalized overlapping integrals between FM and HM. Using 
known expressions for the transverse modal functions of the step-index fiber || , we have 
calculated them numerically In Fig. la, we display 70 and 71 vs. the waveguide parameter 

Note that Eqs. (Q) - (||) do not contain four-wave mixing (FWM) terms. Some of 
them might originate from the terms ~ Uq (u±j and its complex conjugate in the model's 
Hamiltonian density. However, the full expressions to be inserted into the Hamiltonian 
are multiplied by the modal angular dependences exp(=|=2«0), hence they will give zero 
upon the angular integration. Another possible source of FWM terms in Eqs. (||]) - (HD 
could be the term ~ UqU* + u*_ and its complex conjugate in the Hamiltonian density. In 
these expressions, the angular dependence cancels out, hence the angular integration will 
not nullify them. However, the corresponding terms in Eqs. @ - (|5]) will be rapidly 
oscillating in z because of the difference between the propagation constants k and k\. A 
straightforward consideration yields an estimate for the relative wavenumber mismatch 
between FM and HM in the region of interest, \k\ — k \/k « 0.4 (n co — n c \). Taking the 
same estimate for the refractive index difference, n co — n c \ ~ 0.01, as above, we conclude 
that \ki — ko\/k ~ 0.005, which corresponds to the beat length ~ 200 wavelengths. As it 
is disparately small in comparison with any propagation distance relevant to the solitons, 
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all the FWM terms may be neglected. 

As for the dispersion coefficients in Eqs. (1) - (3), their parts accounted for by the 
waveguide geometry can also be calculated for the two modes on the basis of the data 
available from the linear-propagation theory ||. The result of the calculation is shown 
in Fig. lb. It is noteworthy that the waveguide-geometry part of (3q changes its sign. 
One should, however, keep in mind that the full dispersion also contains a material (bulk) 
contribution, that may be essentially larger than that displayed in Fig. lb. 

Thus, the analysis of the interaction between solitons must admit different (but both 
negative, i.e., anomalous 0) effective dispersions (3q and (3\ in Eqs. © - ©• Together 
with the nonstandard values of the XPM coefficients 70 and 71, these features constitute 
an essential mathematical difference of the present model from a usual three-channel 
WDM one (see, e.g., Ref. 0]). 

The fundamental and helical modes are also characterized by a difference in their 
group velocities, which plays a crucially important role in the analysis of soliton-soliton 
collisions. Continuing the above estimates of the physical parameters for the present case, 
we obtain 

(^r) m odeAv~5-l(T 6 (6) 

for the relative group-velocity difference between the modes. As concerns the possibility 
to use the mode-division doubling of the channels in the WDM system, it is relevant to 
mention that, in the WDM system implemented in the standard telecommunications fiber 
with the dispersion (3 ~ —20 ps 2 /km at 1.54 /im, the relative group- velocity mismatch 
between the adjacent channels is estimated to be 

/v gI ~ 1(T 2 ■ (SX/X) , (7) 
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SX being the wavelength separation between the channels. The case of practical interest 
is AA ~ 1 nm, hence we conclude that the relative group-velocity differences (Q) and (0) 
are on the same order of magnitude. On the other hand, in the dispersion-shifted (DS) or 
dispersion-compensated (DC) fibers, the effective value of the dispersion is much smaller 
than the above-mentioned value —20 ps 2 /km, hence in these cases the inference is that 
the corresponding WDM relative difference is negligible as compared to that between the 
fundamental and helical modes, 

( 5 %)S C) / V P < (^gr)mode ( 8 ) 

For typical solitons to be used in telecommunications, with the temporal width T ~ 10 
ps, the above estimate (5f gr ) mode /v gr ~ 5 ■ 10~ 6 implies that a collision between the FM 
and HM solitons takes place at a propagation distance z co ii ~ T (^gr) mo( je A'gr ~ 500 
m, which is much shorter than any soliton's length scale. This circumstance allows us 
to treat the XPM-mediated interaction as a small perturbation in the course of the fast 
passage of one soliton through the other, as it was done in other contexts in Refs. j7], [S| . 
Note that in the laboratory experiments with subpicosecond solitons, the collision length 
may be < 50 m, implying that the experimental study of the collisions should be quite 
possible. 

However, there is no group- velocity difference between the two HMs u±, hence the 
collision distance for the corresponding solitons may be very large, giving rise to a strong 
crosstalk between them. Moreover, the collision between two solitons with the helicities 
S — ±1 may result in their annihilation or transformation into a pair of S = solitons, 
while, due to the conservation of the topological invariant, the collision between the 
solitons with S = and S = 1 is expected to be much closer to an elastic one. In view of 
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this, it makes sense to assume only the doubling of the number of channels by means of 
the "mode-division multiplexing" (i.e., to use only one HM) in the context of the WDM 
systems, but not tripling, that might seem possible due to the existence of two HMs with 
S = ±1. Irrespective of that, a study of collisions between the solitons with S = +1 and 
S = — 1 is a challenge for experiments with narrow solitons in optical fibers. 



3 Analytical Treatment of Soliton Collisions 

Proceeding to the perturbative analysis of the collision between the solitons carried by FM 
and HM, we should take into regard that, in view of the asymmetry between Eqs. (|3D and 
(|4j), (|5p, the FM and HM solitons may have different widths, To and T\. This circumstance 
makes it technically impossible to base the perturbative treatment of the collision on the 
exact unperturbed soliton waveforms of the sech type, as the corresponding overlapping 
integrals will be intractable. The only possibility to develop an efficient perturbation 
theory is to use, as the zero-order approximation, the Gaussian ansatz for the unperturbed 
soliton solutions to the uncoupled equations (0) and (§), 



u? , \z,T) = A l exp( i K l z-^^y ^ = c; J = 0, 1 , (9) 



where a relation between the amplitude and width of the soliton can be found by means 



of the variational approximation [12 



4 a = V2|#|/7? (10) 

(the propagation constants Ki will not be needed here). In fact, a difference between the 
approximate soliton's shape given by Eq. @ and the exact sech shape is fairly small, see, 
e.g., Fig. 5 in Ref. |12| . 
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A soliton moving in the given reference frame is obtained from Eq. as its Galilean 
transform, 

u t (z, r) = u\°\z, t - ti(z)) exp (-iu t T + iqz) , (11) 

where u>i is an arbitrary transform-generating frequency shift, the propagation-constant 
shift q is not essential, and (cf. Eq. (g)) 

^ = Q-|Ah- (12) 

az 

If the XPM term in Eqs. (1) - (3) is, effectively, a small perturbation (in the case of 
a fast collision, see above), the collision between the solitons may be described as that 
between two quasiparticles interacting through an effective potential. Following the lines 
of the analysis developed for similar problems earlier [0, § , it is straightforward to derive 
the following perturbation-induced evolution equations for the solitons' frequency shifts: 



dui 4|$|7z d 



i H / ' 17 ex P 

d z T^JTg + T? Mi 



(f i - 1, 



o 



\2 



0, (13) 



2(1? + 2?) 

where Eq. (|TD|) was used to eliminate the amplitudes in favor of the widths T\ (recall that 
7z are the relative XPM coupling constants in Eqs. (Q) - (|5])). Combining Eqs. (|I3|) with 
Eqs. (|i~2l) and assuming, in the first approximation, T\ = const furnish a closed dynamical 
system governing the evolution of the temporal positions ti of the two solitons. 

To further apply the perturbation theory to Eq. ([13]), we recall that, according to 
the estimates obtained above, the difference of the inverse group velocities, c = c\ — cq, 
is, effectively, a large parameter. Hence, in the lowest-order approximation, one may set 
t 1 —t m cz in the argument of the exponential in Eq. (0), thus strongly simplifying the 
equation: 



dun _ 4(-l) / |A| 7z d <r; 



dz cT^Tg + T? dz 
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cxp 



2 



2(t 2 + n 



(14) 



To proceed, it is necessary to specify a type of the collision to be considered. One 
should distinguish between "complete" and "incomplete" collisions . In the former case, 
the solitons are, originally, far separated; in the course of the interaction, the faster soliton 
catches up with the slower one and passes it. In the first approximation, the complete 
interaction does not give rise to a net frequency shift (a change of the frequency would be 
tantamount to a change of the soliton's velocity, according to Eq. fll2|)), as the integration 
of the right-hand side of Eq. (14) from z = — oo to z = +00 yields zero. However, finding 



a nonzero instantaneous frequency shift from Eq. fllil) , inserting it into Eq. fljjD, and 
integrating the latter equation yield a nonzero collision-induced position shift Sti of the 
soliton's center, which is the main effect of the complete collision. A final result can 
be conveniently written as a relative position shift, normalized to the soliton's temporal 
width: 

* (-ir-ws-gU (i5) 

Ti c 2 T Ti 
An "incomplete" collision takes place if the solitons are essentially overlapped at the 
initial point, z = 0. This kind of the collision is more significant, as it gives rise to 
a nonzero net frequency shift 5ui (hence, to a velocity shift too). The most important 
(dangerous) case is that when centers of the colliding solitons exactly coincide at z — 0. In 
this case, dui is found by straightforward integration of Eq. ( |T4"D from z = to z = +00. 
The result can be presented in a more natural form, multiplying the net frequency shift by 
the soliton's temporal width (i.e., normalizing the frequency shift to the soliton's spectral 
width): 

Tfa = 4(-l)^ m ' 7 , . (16) 

cT^y/Tg + Tg 

The only difference of Eqs. ( [To] ) and ( |T6"D from similar results for the usual WDM 
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system are the specific XPM coefficients 7;, which are = 1 in the usual case. The most 
promising range for the applications is around V = 3.6 (Fig. la), which gives 

70 « 0.98, 7! « 0.62. (17) 

This implies that the crosstalk between the FM and HM solitons is attenuated by the 
factor 0.62 for the HM soliton, as compared to the usual WDM system, while for the 
FM soliton the crosstalk strength is not different from that in the usual system (taking, 
instead, the values around V = 2.4 (Fig. la), we will get small 71 ~ 0.25, but large 
7o ~ 1.66). 

Thus, if the set of the HM and DM modes is used to double the number of the 
channels in the WDM system, we conclude that the crosstalk for the HM solitons due to 
their collisions with the FM ones in any channel is attenuated, against the usual strength, 
by the above-mentioned factor « 0.62. 

Note that the conclusion concerning the comparison with the WDM crosstalk pertains 
to the case when the WDM system is realized in the standard telecommunications fiber: as 
it was concluded above, in this case the group- velocity difference between the fundamental 
and helical modes is on the same order of magnitude as the group- velocity mismatch 
between adjacent WDM channels, see Eqs. (||) and (0). On the contrary to this, in 
the DS or DC fiber link the group-velocity difference between the FM and HM channels 
is much larger than that between the WDM ones, hence the FM-HM crosstalk is much 
weaker than between the WDM channels, according to Eqs. fll5]) and (IB). 
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4 Numerical Simulations of the Collision 

Direct simulations of the soliton collisions within the framework of Eqs. (|3]) and (|J) for the 
two modes demonstrate that, although the above-mentioned frequency-shift-attenuation 
factor 0.62 is not really small, sometimes it may be important. In Fig. 2a, we display an 
example of a disastrous incomplete collision in the usual WDM system, which leads to a 
merger of the solitons into a breather (the simulations of the complete collision at the same 
values of the parameters shows that it is fairly mild, producing only small position shifts 
of the solitons). Replacing the usual XPM coefficients 7 = 1 by those for the FM-HM 
collision, given by Eq. (|17]), we find that the same solitons survive the incomplete collision 
(Fig. 2b). 

5 Conclusion 

We have considered solitons in a nonlinear optical fiber in a parametric region where the 
fiber supports exactly two distinct modes, the fundamental one and the first helical mode, 
which allows one to double the number of soliton channels in the fiber. We have intro- 
duced a system of coupled NLS equations for the two modes and computed nonstandard 
values of the relative XPM coupling constants in it. Then, we investigated analytically 
and numerically results of both "complete" and "incomplete" collisions between solitons 
carried by the helical and fundamental modes, concluding that the crosstalk effects are 
partly attenuated for the helical solitons. The crosstalk between the two modes is found to 
be additionally suppressed in comparison with that in the WDM system in the case when 
the latter system is realized in a dispersion-shifted or dispersion-compensated fiber. Irre- 
spective of the possible applications, experimental observation of narrow helical (spinning) 
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temporal solitons in ordinary or hollow optical fibers is a challenge, especially because it 
has been recently demonstrated that spinning spatiotemporal solitons ("light bullets") 
are unstable in bulk media. 
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FIGURE CAPTIONS 

Fig. 1. The normalized XPM coefficients (a) and the normalized waveguide dispersion 
coefficients (b), vs. the waveguide parameter V for the fundamental (S = 0) and he- 
lical (S — 1) modes in the standard step-index fiber. The full waveguide dispersion is 
(n c \VA/Xc) D, where c is the light velocity, and A = (n co — n c \)/n co . In the plot (a), the 
dashed line shows the usual value (= 1) of the XPM coefficient. 

Fig. 2. Collisions of initially overlapped solitons with the amplitudes |w| ma x = Mmax = 
2 and the relative inverse group velocity c\ — c = 29 (recall it must be a large parameter) 
simulated with Eqs. (1) and (2), in which f3 Q — fi x — — 1: (a) the usual WDM system, 
with 7 = 7x = 1; (b) the fundamental and helical solitons, with 7 = 0.98 and 71 = 0.62. 
The simulations were carried out from z = up to a distance equal to six soliton periods, 
which is z = 6. 
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Fig. 1(b) 




Fig.2(b) 



